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ABSTRACT 

We present the faint end of number counts at 1.3 mm (238 GHz) obtained with the Atacama Large 
Mihimeter/submihimeter Array (ALMA). Band 6 observations were carried out targeting 20 star- 
forming galaxies at z ~ 1.4 in the Subaru/XMM-iYewton Deep Survey field. In the observations, we 
serendipitously detect 15 sources (>3.8cr, 5'i.3 mm = 0.15-0.61 mJy) other than the targeted sources. 
We create number counts by using these 'sub-mJy sources', which probe the faintest fiux range among 
surveys at millimeter wavelengths. The number counts are consistent with (fiux-scaled) number counts 
at 850 /im and 870 /im obtained with gravitational lensing clusters. The ALMA number counts agree 
well with model predictions, which suggest that these sub-mJy populations are more like 'normal' 
star-forming galaxies than 'classical' SMGs with intense star-forming activity. In this fiux range, 
~80% of the extragalactic background light at 1.3 mm is resolved into individual sources. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — galaxies: 
ISM — cosmology: observations — submillimeter: galaxies 



1. INTRODUCTION 

Extragalactic surveys at millimeter and submillime- 
ter have uncovered a population of dusty star-forming 
galaxies at high redshifts, so-called submillimeter galax- 
ies (SMGs). SMGs are dusty starburst galaxies with 
star-formation ra tes (SFRs) of 10^-10^ Mq yr"^ (e.g., 
iBlain et 'al]l2002f ). SMGs are important for understand- 
ing of the process of galaxy formation, the cosmic star- 
formation history, and the origin of extragalactic back- 
ground light (EBL). The EBL is the integral of unre- 
solved emission from extragalactic sources and contains 
information on the history of galaxy evolution. A large 
fraction of the EBL at millimeter and submillimeter 
wavelengths is thought to be contributed by high-redshift 
dusty galaxies (Lagache et al.l l2005f ). Previous blank 
field surveys reso l ved 20%-40% of the EBL at 850 /im 
(e.g., lEales et al.l 120001: iCoppin et al 1 12006|) and __10%- 
20% at 1 mm (e.g .,'Greve et al."2004'; 'Perera et^|2QQ8|; 
iScott et al.ll2010r Hatsukade et al. 2011). Meanwhile, in 
cluster fields 50%-100% of the EBL are res olved by uti- 
lizing the effect of yavitational lensing (e ^., lSmail et al.l 
119971 : iCowie et al.ll2002t iKnudsen et al.l 12008). suggest- 
ing that a large fraction of the EBL originates from faint 
sources {Si mm < 1 mJy) of which number counts are not 
yet well constrained due to the limited sensitivity and 
confusion limit for the single-dish surveys. Compared 
to SMGs detected in previous surveys, fainter submil- 
limeter sources are possibly connected to more general 
star- forming galaxies, or galaxies with hotter dust. Ob- 
servations of fainter sources are needed to understand the 
formation and evolution of galaxies. 

The Atacama Large Millimeter/submillimeter Array 
(ALMA) enables us to explore much fainter sources with- 
out the effect of confusion limit because of its high sen- 
sitivity and high angular resolution. In this letter, we 
present the faint end of number counts at 1.3 mm re- 
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vealed by ALMA band 6 observations. The arrangement 
of this paper is as follows. Section [2] outlines the obser- 
vations and data reduction. In Section [3l we describe 
the method of source extraction and simulations carried 
out to create number counts. In Section HI we compare 
the number counts with other observational results and 
model predictions, and estimate the contribution of the 
resolved 1.3 mm sources to the EBL. A summary is pre- 
sented in Section [5l 

2. OBSERVATIONS AND DATA REDUCTION 

We conducted ALMA band 6 observations to- 
ward 20 star- forming galaxies at z ~ 1.4 in the 
Subaru XM M/Newton Deep Survey (SXDS) field 
(|Furusawa et al.ll200 8). The targets were extracted from 
a stellar mass limit sample whose redshifts were ob- 
tained by near-infrared (NIR) spectroscop y with the Fi- 
bre Multi-Object Spectro-grap h (FMOS; iK imura e t al.l 
I2010f ) on the Subaru telescope ( lYabe et al.l [2012. 2011). 
The observations were carried out in August 9th, 11th, 
15th, and 26th, 2012 with 23-25 antennas in the ex- 
tended configuration during the ALMA Cycle session. 
The correlator was used in the frequency domain mode 
with a bandwidth of 1875 MHz (488.28 kHz x 3840 
channels). The targets were divided into four groups 
with different spectral settings to observe the redshifted 
CO (5-4) line as efficiently as possible. The local oscil- 
lator frequencies for the four groups were 231.198 GHz, 
236.168 GHz, 240.380 GHz, and 244.166 GHz (typical 
frequency of 238 GHz). Three or four spectral windows 
with a bandwidth of 1875 MHz, each with a different 
central frequency, were used in each group to cover the 
CO (5-4) line and simultaneously detect the dust thermal 
continuum emission. Ganymede, Callisto, and Uranus 
were observed as a flux calibrator. The bandpass and 
phase were calibrated with J2253+161 and J0204— 170, 
respectively. We obtained 20 pointings centered on the 
20 targets, each with on-source observing time of 8-15 
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Fig. 1. — Average number of positive and negative peaks in a 
map as a function of SN threshold. 

minutes. The full width at half maximum (FWHM) of 
the primary beam is 26'^ 

The data were reduced with the Common Astronomy 
Software Applications (CASA; Mc Mullin et~all \2007h 
package in a standard manner. The maps were pro- 
cessed with the CLEAN algorithm with the natural weight- 
ing, which gives the final synthesized beamsize of ^0.6"- 
1.3''. The rms noise level of the 20 maps is 0.04-0.10 
mJy beam~^ and is almost constant in each map un- 
corrected for the primary beam attenuation. In what 
follows, we use the area within the primary beam of the 
20 maps. In this letter, we focus on serendipitously de- 
tected sources other than the targeted sources to exam- 
ine the source number counts. Details on the targeted 
star- forming galaxies are discussed elsewhere. 

3. SOURCE EXTRACTION AND NUMBER COUNTS 
3.1. Source and Spurious Detection 

Source extraction is conducted on the 20 maps (before 
primary beam correction) where all the sources with a 
peak signal-to-noise ratio (SN) above 3.5 are CLEANed. 
We adopt the detection threshold of 3.8cr, which is de- 
termined from the number of spurious sources expected 
in the maps. The number of spurious source is estimated 
by counting the negative peaks in each map. Spurious 
source extraction is performed on the maps multiplied 
by —1. The spatial distribution of the spurious sources 
is almost uniform in each map. Fig. [1] shows the number 
of positive and negative peaks in a map as a function of 
SN threshold averaged over the 20 fields. We find that 
the number of spurious sources is less than one at SN 
>3.8. 

We detect 15 sources (excluding the target sources) 
in the 20 maps with the raw flux density (uncorrected 
for the primary beam attenuation) of 0.15-0.53 mJy 
(3.86-5.84cr). All the sources are point sources or only 
marginally resolved. A few sources are possibly detected 
in emission lines which may dominate their flux density. 
Examples of detected sources in the maps with different 
SNs are shown in Fig. [2] (top). Fig. [2] (bottom) shows the 
source positions in the maps. No clear clustering feature 
is seen in the sky distribution of the sources. 
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Fig. 2. — Top: Examples of detections in four maps (5^^ x 5^^). 
Contours start from ±3. Oct with a 0.5cr step. Negative contours are 
presented in dashed lines. The synthesized beam is shown in the 
bottom left in each map. Bottom: Positions of the 15 source in the 
maps. The size of each symbol is scaled by signal-to-noise ratios. 
Dashed circle represents the primary beam. 



3.2. Completeness 

We calculate the completeness, which is the rate at 
which a source is expected to be detected in a map, to 
see the effect of noise fluctuation on the source detection. 
Because the source extraction is conducted in the maps 
uncorrected for primary beam attenuation, the complete- 
ness calculation for each source is also conducted in the 
same map by using the raw flux density. We subtract 
>3.5cr sources from each map, and inject a flux-scaled 
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Fig. 3. — Effective area as a function of flux density, where a 
source with a flux density (corrected for primary beam attenuation) 
can be detected with >3.8cr. 
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synthesized beam into the map. Because the rms noise 
level is almost constant in each map, the input position 
is randomly selected in the map. When the input source 
is extracted within l.O'Oof its input position with >3.8cr, 
the source is considered to be recovered. We repeat 
this 1000 times for the raw flux density of each source. 
The completeness for the 15 sources ranges from 77% to 
100%. 

In the course of this simulation, we calculate the ratio 
between the input and output fluxes to estimate the in- 
trinsic flux density of the detected sources. The flux ratio 
calculated for each source in this simulation is 1.00-1.14. 
The flux density of the sources corrected for this effect 
and the the primary beam attenuation is 0.15-0.61 mJy. 

3.3. Effective Area 

Because the rms noise levels are different in the 20 
maps and the primary beam response in a map changes 
depending on the distance from the map center, the sur- 
vey area for each source should be calculated depending 
on its intrinsic flux density (corrected for primary beam 
attenuation). We calculate the effective area against the 
intrinsic flux density in each map; By using the rms noise 
levels corrected for primary beam attenuation, we find a 
distance from the map center within which the source is 
detected at >3.8<j in the raw map uncorrected for pri- 
mary beam attenuation. 

3.4. Number Counts 

By using the 15 serendipitously-detected sources, we 
create cumulative number counts. To create number 
counts, we correct for the contamination of spurious 
sources, the incompleteness, and the survey area for each 
source. The contamination of spurious sources to each 
source (Cgpurious) is estimated from the average spu- 
rious detections at its SN (Fig. [1]) and is subtracted 
g). Then it is divided by the com- 
ious) /completeness), and finally cor- 
rected for the survey area by using the effective area 

^ Note that the effect of the choice of 0.5^^ or 2^^ on the number 
counts is within the errors. 



from 1 (1 — Cs] 
pleteness ((1 — Cg] 



Fig. 4. — Cumulative number counts at 1.3 mm obtained in 
the ALMA band 6 observations (red squares for >3.8cr sources 
and blue circles for >4.0cr sources). For com parison, we plot 
the ALMA 870 /xm counts (IKarim et al.l 120121 ) in the Extended 
Chandra Deep Field South, A zTEC 1.1 m m combined counts 
of 6 blank-flelds surv eys (Sco tt et al.|[2QTl ). LABOCA 870 //m 
(Johansson et al. 2011) and SCUBA 850 /xm counts around galaxy 
clusters (Knudsen et al. 200 ^ . We also show the model predictions 
of Shi mizu et al.l (l2012l). iHavward et al.l (|20T3h . IBethermin et al.l 
(2012). and Ida Cunha et al.l (|2013l ). The 1.1 mm, 850 fim, and 
870 /im counts are offset to an equivalent 1.3mm flux density using 

SCahng factors of S'l.Smm/'S'l.lmm = 0.71, S'l.3mm/*S'870/^m = 0.38, 
and S'l.Smm/'S'sSO/xm = 0.36. 



TABLE 1 

Number counts at 1.3 mm derived from the >3.8a sources. 



s 

(mJy) 
(1) 


A/source 

(2) 


N(> S) 

(xl03 deg-2) 

(3) 


0.15 

0.24 
0.37 
0.59 


2 

4 
7 
2 


ro 1+108 
18-3+691 

091+13:0 

^•^1-46. 6 



Note. — The errors are Poisson confldence limits of 84.13%, 
corresponding to la for Gaussian statistics. (1) Flux density; (2) 
Number of sources in the bin; (3) Cumulative number counts. 



((1 — Cspurious) /completeness/effective area). The up- 
per and lower errors are calculated for each source from 
Poisson confidence limits of 8 4.13%, which correspond 
to la for Gaussian statistics (|Gehrelslll986f ). and then 
corrected for the incompleteness and survey area in the 
same manner as described above. The derived number 
counts are shown in Fig. |4] and Table [TJ We also plot 
the number counts created by using 6 >4.0<j sources to 
see the effect of the detection threshold. Both counts are 
consistent within the margin of errors. 

4. DISCUSSION 

By using the "sub-mJy sources" detected with ALMA, 
we derived cumulative number counts, which probe cur- 
rently the faintest end at 1 mm wavelengths. We com- 
pare the number counts with other observational results 
and model predictions. 
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4.1. Comparisons with Previous Results 

We compare our number counts with previous ob- 
servational results. In Fig. IH we plot number counts 
obtained by SCUBA 850 um (Knudsen et al._ J008), 
LABOC A 870 um (Ijohansson et all I2QTT1). AzTEC 
1.1 mm CScott et al." 2012D and A LMA band 7 (870 /im) 
observations (Karim et al.l [20121 ). The flux density of 
these counts are scaled to 1.3 mm flux density by us- 
ing a modifled black body with typical values for SMGs; 
spect ral index of /3 = 1.5 and dust tem pera ture of 35 K 
(e.g., iKovacs et atl 120061: ICoppin et al.l 12008.) locat ed at 
z = 2.5 fe.g. JChapman et al.ll2005l : lYun et al]l2012t ). 

The faintest counts are probed by using gravitationally 
lensing clusters at 850 /im (Knudsen et al. 2008) and at 
870 /im (j Johansson et al.ll201l[ ). The ALMA 1.3 mm 
counts are consistent with those counts, suggesting that 
b oth wavelengths tr ace the same population at faint end. 

iScott et arr(|2012[ ) derived currently the most reliable 
number counts at 1.1 mm by combining data in six 
blank flelds surveyed by AzTEC. The 1.3 mm counts 
appear to be consistent with those at 1.1 mm, except 
for the faintest bin of the 1.1 mm counts. They sug- 
gest a turnover in the Euclidean-normalized counts at 
*5'i.i mm ^ 2 mJy. However, our data points at fainter 
flux range do not show such turnover. We note that 
number counts derived from targeted observations may 
show an exce ss of source c ou nts. 

Recently, iKarim et al! (J2012D conducted ALMA 
870 /im observations of 122 sources detected in the 
LABOCA Extended Chandra Deep Field South Submil- 
limeter survey. The high-resolution observations pin- 
point the source positions, and some of the LABOCA 
(FWHM ~ 19'') sources are resolved into multiple 
sources. Number counts are derived from the ALMA 
results without the effect of source confusion. Although 
the extrapolation of the fainter end of the 870 /im num- 
ber counts appears to be lower than the 1.3 mm counts, 
deeper observations at the same depth in flux density 
as the flux ranges of 1.3 mm observations are needed to 
compare both results quantitatively. 

4.2. Comparisons with Model Predictions 

We compare the 1.3 mm counts with recent model 
predic tions with different appro aches b y iShimizu et al. l 
(120121) . iBethermin et al.' (2012), and "da Cunha et alj 
(|2013f ). Fig. m illustrates that the three models agree 
w ell with the ALMA 1 .3 mm number counts. 

IShimizu et al.l (|2012l ) perform large cosmological hy- 
dro dynamic simulations and simulate the properties of 
SMGs by calculating the reprocessing of stellar light by 
dust grains into far-IR to millimeter wavelengths in a 
self-consistent manner. Their model reproduces number 
counts and clustering amplitude of SMGs. The ALMA 
number counts are consistent with the predicted number 
counts calculated for 238 GHz with this model. In the 
flux range of 0.1-0.6 mJy, the simulated galaxies have a 
stellar mass of -5 x 10^-10^^ M© and an SFR of -10- 
100 Mq yr-^ 

The model of IBethermin et all (|2012[ ) is based on the 
redshift evolution of the mass function of star-forming 
galaxies, speciflc SFR distribution at flxed stellar mass, 
and spectral energy distributions (SEDs) for the two star- 
formation modes (main-sequence and starburst). The 



main-sequence galaxies, which follow a tight correlation 
between stellar mass and SFR, have steady star forma- 
tion, while starburst galaxies, which are outliers to the 
correlation, h ave hi^h specific SF Rs possibly driven by 
mergers (e.g., iNoeske et all 120071 ). Their model repro- 
duces IR to radio number counts and counts per red- 
shift slice of Herschel observations. At fainter flux den- 
sity (5'i.imm ^ 1 mJy) main-sequence galaxies dominate 
the number counts by an order of magnitude, suggesting 
that the ALMA sub-mJy sources are more like normal 
star-forming galaxies than classical SMGs with intense 
st ar- forming activ i ty. 

I da Cunha et al.l (|2013 i) create self-consistent models of 
the observed optical/NIR SEDs of galaxies detected in 
the Hubble Ultra Deep Field. They combine the attenu- 
ated stellar spectra with a library of IR emission models, 
which are consistent with the observed optical/NIR emis- 
sion in terms of energy balance, and estimate the contin- 
uum flux at millimeter and submillimeter. The predicted 
number counts agree with the previous observed number 
counts and extragalactic background light. They simu- 
late deep ALMA observations with a continuum sensi- 
tivity of IcTcontinuum = 5.1 X 10~^ mJy. In their simu- 
lation, the redshift distribution of detected galaxies in 
the ALMA band 6 is relatively flat up to z = 5. The 
detected galaxies typically have a stellar mass of 10^- 
10^^ Mo and an SFR of a few x (0.1-10) M© yr-\ which 
are more massive, more star-forming, and more dusty 
galaxies among their optical/NIR-detected galaxy sam- 
ple. The agreement between the model and our ALMA 
number counts implies that the sub-mJy sources are mas- 
sive and high-SFR end of optical/NIR-detected galaxies. 

4.3. Contribution to the Extragalactic Background Light 

We calculate the fraction of the EBL revolved into dis- 
crete sources in the ALMA 1.3 mm observations. The to- 
tal flux density of the >3.8cr sources, which is applied all 
the corrections conducted in the derivation of the number 
counts (Section[3|), is ~13 Jy deg~^. The EBL at 1.3 mm 
based on the measurements by the C osmic Background 
Explorer satellite is 14-18 Jy deg-^ (jPuget et all 119961 : 
iFixsenet al1ll998D . Therefore, we have resolved -80% 
of the EBL into discrete sources, which is signiflcantly 
higher compared to results in the previo us deep sur- 
veys (10%-20%) at 1 mm wavelengths (e.g.,lGreve et"all 
120041 : IScott et al. 2010; Ha tsukade et al.l 120 1 11 thanks 
to the high sensitivity and high angular resolution of 
ALMA. The index of a power-law flt to the number 
counts is steeper than —1, indicating that the contri- 
bution of fainter sources is signiflcant. 

5. SUMMARY 

We conducted deep 1.3 mm observations with ALMA 
toward 20 star- forming galaxies at 2: — 1.4 in the SXDS 
fleld. We serendipitously detected 15 sources (>3.8cr) 
other than the targeted star- forming galaxies. The flux 
densities of the ALMA sources are 81,3 mm = 0.15- 
0.61 mJy, which are an order of magnitude fainter than 
those of SMGs detected in previous surveys. By us- 
ing these 'sub-mJy sources', we created number counts, 
which probe the faintest flux range at millimeter wave- 
lengths. The number counts are consistent with flux- 
scaled 850 /im and 870 /im counts obtained with gravita- 
tional lensing clusters, which trace the faint end of num- 
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ber counts at submillimeter wavelengths. The number 
counts agree well with model predictions, which suggest 
that the sub-mJy sources are moderately star-forming 
galaxies with SFR ^100 Mq yr~^ rather than classi- 
cal SMGs with intense star formation. We find that 
our ALMA observations resolves ~80% of the EBL at 
1.3 mm. 
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